Variational perturbation approach to the Coulomb electron gas 
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The efficiency of the variational perturbation theory [Phys. Rev. C 62, 045503 (2000)] formulated 
recently for many-particle systems is examined by calculating the ground state correlation energy of 
the 3D electron gas with the Coulomb interaction. The perturbation beyond a variational result can 
be carried out systematically by the modified Wick's theorem which defines a contraction rule about 
the renormalized perturbation. Utilizing the theorem, variational ring diagrams of the electron gas 
are summed up. As a result, the correlation energy is found to be much closer to the result of the 
Green's function Monte Carlo calculation than that of the conventional ring approximation is. 



I. INTRODUCTION 

A variational method has been widely used as a con- 
venient and powerful tool to calculate physical quantities 
such as energy and order parameters [1-5] . The most im- 
portant advantage of this approach consists in its appli- 
cability even for the strongly correlated systems through 
suitable trial wavefunctions. However, it is extremely dif- 
ficult to improve the results for the systems with many 
degrees of freedom. In contrast, results of the perturba- 
tive approach are systematically improved by calculating 
the further corrections using the Feynmann rule [6-8]. 
However, the convergency can not be guaranteed in the 
unweakly coupled cases. 

There have appeared many efforts trying to combine 
the above two approaches in order to overcome each 
drawback [9-14]. Many studies have been performed 
mainly in the areas of relativistic field theories and quan- 
tum mechanics under the principle of minimal sensitivity 
(PMS) [9], which says that approximated physical quan- 
tities in the perturbation theory should be minimized 
with respect to the parameters absent in the original 
Hamiltonian. These studies are called optimized pertur- 
bation theory or variational perturbation theory (VPT), 
which is applied very successfully to quantum mechani- 
cal problems such as the anharmonic oscillator and the 
double- well potential. However, we are faced with ex- 
tremely complex higher-order calculations in the case of 
quantum field theories. 

Recently, papers on another kind of VPT have been 
published [15-18]. This approach is not related to the 
PMS, because the minimization is carried out at the ze- 
roth order and the variational parameters are fixed before 
the perturbative calculation. It has a simple expansion 
rule, which makes the VPT approach to many-particle 
systems more manageable. In this paper, the VPT for 
fcrmion systems of ref. [17] is briefly reviewed, and as 
a test on the efficiency of the method, we calculate the 
ground state correlation energy of the three dimensional 
Coulomb electron gas within a ring approximation. 

II. VARIATIONAL PERTURBATION THEORY 

In this section, we briefly review the VPT in the func- 
tional integral formalism [17]. The partition function for 



fermion systems with two-body interaction is generally 
given by 



Z = / D[i/nl>] exp(- 
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where tpi ,ipi are Grassmann variables and each subscript 
stands for all possible quantum numbers including an 
imaginary time. The integration and summation symbols 
are omitted by the summation convention. Here, G° is 
the bare Green's function matrix and V is an interaction 
tensor. The second term of the exponent in Eq.(l) is 
the perturbative term in the conventional perturbation 
theory. Now, introducing a variational Green's function 
G, the exponent is rewritten as 

- Gr>^. - (G% x - Gt. 1 )^ - Vi m $iMk1>i ■ (2) 

Then, the second and third terms in Eq.(2) are regarded 
as a renormalized interaction. According to the Jensen- 
Peierls' inequality [19,20], we get the inequality between 
the thermodynamic potential Q and the variational ther- 
modynamic potential Q(G); fl < fi(G). Minimizing 
f2(G) with respect to G, we obtain an equation of G as 



G - G° 
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where (6|S|a 

corresponds to the self-consistent Hartree-Fock Green's 
function for the interacting fermion system. Using the 
condition Eq.(3), the minimized variational thermody- 
namic potential n m i n and the correction part Ail = 
— iljnjn are arranged as follows, according to the nota- 
tion in ref. [17]; 
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(4) 
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G 



represents the thermal 



average using G and the subscript G indicates the con- 
nected contractions among all possible contractions by 
the modified Wick's theorem Eq.(7) below. During the 
derivation of the above equations, Grassmann variables fj, 



rj are introduced as source fields of the original fields ip, ip, 
and the renormalizcd interaction of Eq.(2) is expressed as 
functional derivatives about source fields. After integrat- 
ing out the Gaussian integral (/ D[^pip] cxp(— -i/JiG" 1 ^ + 
fji'ipi+'ipirji) — DetG~ 1 e'' iG ^' ?J ') and rearranging the ther- 
modynamic potential under the condition Eq.(3), we ob- 
tain Eqs.(4) and (5). Here, the primed derivative opera- 
tion (• • •)' is defined as 
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Carrying out the nth order calculation of the above op- 
eration, we find modified Wick's theorem for the pertur- 
bative expansion of the present VPT as 



: Sum of all possible contractions 
between different cells. 



0< 



(7) 

Here we call one (• • •)' a unit cell and the contraction 
rule is as follows; each d Vi is paired with a d flj in a dif- 
ferent cell and then, each pair of derivatives d m dfjj is 
replaced by Green's function Gji. In this pairing, every 
permutation of two derivatives changes the sign. The 
differences from the original Wick's theorem are to for- 
bid the intracell contraction and to use the renormalized 
Green's function G. The unit cell with an interaction, 
Vijki (dmdrjjd^dfj,) , can be described by a diagram as 
in Fig.l (a). A contraction is to join two incomming 
and outgoing lines together and to assign it a propaga- 
tor G which corresponds to a double line. In Fig.l (b), 
we show ring diagrams with renormalized propagators, 
where the first order diagram does not exist because it is 
an intracell-joining diagram. 

For example, the second order perturbation of Aft is 
given by; 
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III. COULOMB ELECTRON GAS 

The present formulation can be easily applied to any 
fermion systems with two-body interaction. In this sec- 
tion, the ground state correlation energy of the three di- 
mensional Coulomb electron gas is calculated using the 
variational ring approximation which is described by ring 
diagrams in Fig.l (b). The model Hamiltonian is given 

by 

H = Y^ ^a\ a a ka + i ' v {<l)aia a i>a> a k>-q<y>ak+ q <y , 



ka 



kk'qaa' 



(a) 






(b) 



FIG. 1: (a) Unit cell to construct Feynmann diagrams; (b) 
the summation of variational ring diagrams, the first order 
diagram cannot exist in the present theory. 



where a\ a (a ka ) is a creation(annihilation) operator of an 
electron with wavevector k and spin a. We have defined 
e k = h 2 k 2 /2m e and v(q) = 47re 2 /Vq 2 . Here, V, m e and 
e are the system volume, the electron mass and charge 
respectively. The primed summation indicates that the 
q = term is excluded because of the cancellation with 
the positively charged background. In this model, the 
interaction Vijki of Eq.(l) is expressed as 

Vijkl > ^ / fcl<TlTi,fc2<T2T2, fc3CT3T3, fc 4 <T4T4 

= -jv{\k± - k 1 \)5 kl+k2yk3+k4 S rTiai 5 a2Cr3 

x 5(n - t 2 )S(t 3 - t 4 )<5(t 2 - r 3 ) . (10) 

The variational Green's function G is determined by 
Eq.(3). For a homogeneous system, G kaT,kar' is writ- 
ten as Gka(r — t') and the Fourier transform of Eq.(3) 
with respect to time gives 



G ka (iuj n ) = - — 
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(11) 



(9) 



where uj n — (2n + l)n/f3 is the fermion Matsubara fre- 
quency with integer n, and e ka = e k - J2 q v (l) n k+qa 
is the renormalized electron energy using the Fermi dis- 
tribution function with the renormalized energy, n ka = 

Therefore, the minimized variational thermodynamic 
potential ^ m i n at T — OK is 

^min = ^(h ~~ ^ Hka + \ ^2 v (<l) n k+qanka ■ (12) 
ka kqa 

This is the Hartree-Fock result. The second term can- 
cels the doubly counted interaction energy in the first 
term to give the singly counted result. If the interac- 
tion part of ha and the second term are summed up, 
f2 m ; n has the same form as the conventional first order 

result fii; VL X = J2ka( e k ~ - 3 T, kq a V (l) n k+qa n ka 

, where n\ a = l/(eP( ek °~^ + 1). In addition, the Fermi 



distribution function with a renormalized energy is a step 
function in the ground state like that with a bare energy. 
Therefore, the Hartree-Fock ground state energy is equal 
to the conventional first order one, which results from the 
spherical symmetry of electron gas [21] . Hence, one might 
think that the ground state energy obtained by higher or- 
der calculations through the present VPT would produce 
the same result as that of the conventional perturbation 
method. However, as we will see below, they are different 
because the VPT propagator has a renormalized energy 
ika and furthermore, there are some forbidden diagrams 
in the VPT expansion, which are allowed in the conven- 
tional perturbation. 

Since the present formalism parallels exactly the con- 
ventional perturbation except for the modification of the 
Wick's theorem and the renormalization of the propaga- 
tor, we can carry out the ring diagram summation with- 
out any difficulty. The ring contribution to the thermo- 
dynamic potential is depicted by Fig.l (b). The double 
line represents a renormalized propagator G. The varia- 
tional ring diagrams are summed up to the logarithmic 
function as in the conventional ring approximation [6], 
namely, 

^ ]T ln[l + v(q)F(q, u n )] - ±- £ v(q)F(q, v n ) , (13) 

where v n = 2irn/ (3 is the boson Matsubara frequency and 
F(q, v n ) is the renormalized Lindhard function, 

F{q,v n ) = Y, - nk p°- nk >Z (14) 

At zero temperature, the discrete frequency summa- 
tion i^p^Yliv ^ S replaced by a continuous integral J dv. 

Figure 2 shows Lindhard functions at zero tempera- 
ture. The small(large) dotted line corresponds to the 
bare(renormalized) Lindhard function. According to the 
linear response theory, the Lindhard function Fo(q,v) 
has important properties that it is proportional to the 
spin and charge susceptibilities of free fermion systems, 
Xo(Qi v ) an d Xoili 1 ')- ln addition, Fo(0,0) is equal to 
the density of states at the Fermi surface, -D(cf) [22]. 
We note that the values of F(q, 0) is much reduced by 
the inclusion of the exchange energy into a bare energy 
band. Furthermore, the behavior near q = is com- 
pletely different, which originate from the much larger 
energy slope around the Fermi surface than the bare one. 
The difference between F(q, v) and Fo(q, v) results in dif- 
ferent ground state energies between the variational and 
conventional ring approximations. 

The chemical potential fi is an independent variable 
of the thermodynamic potential il(/i), hence in order 
to obtain the free energy as a function of electron den- 
sity, we should change the independent variable using 
the Legendere transformation; F(N) — + fiN and 
— §Q — N, where V is the number of particles and F(N) 
is the free energy. However, this relation can not be 
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FIG. 2: The variational Lindhard function F(q, 0) (large dot- 
ted line) at the zero temperature and frequency; the small 
dotted line represents the free electron Lindhard function 
Fo(q, 0). The magnitudes are renormalized by fo(0,0). 

applied directly in approximations except for Hartree- 
Fock's and Baym's self-consistent schemes [23] because 
— |p does not represent the particle number. Instead, we 
use the fact that the Fermi wavevector does not change 
with the inclusion of correlation according to Luttinger's 
theorem [24]. Therefore, we approximate the free en- 
ergy as follows; F(N) ~ ^min(Mo) + Afi rijl g(/xo) + ^ V 
and — dg'^ in | Po = N, where the Fermi wavevector (or 
the density) is fixed by the second relation without the 
correlation part. In the case of the conventional ring ap- 
proximation, if J7 m i n is replaced by a non-interacting part 
fio, the free energy by this transformation reproduces the 
RPA (random phase approximation) correlation energy. 

Figure 3 shows the correlation energy as a function of 
r s which is the averaged relative distance between elec- 
trons defined by V — (4/3)n(a r s ) 3 N, where V, V and 
ao are the system volume, the number of electrons and 
the Bohr radius, respectively. The triangles are the result 
of the Green's function Monte Carlo calculation [25] and 
the better agreement with it is realized by the VPT. The 
conventional rings contain only direct Coulomb interact- 
ing processes without exchange contributions which are 
very important for fermion systems, and thus, the sum- 
mation of them gives, so called, the RPA result [26]. On 
the other hand, the variational rings possess infinitely 
many higher order exchange processes through the vari- 
ational propagator G as well as direct processes. As a 
result, we achieve much improvement in the correlation 
energy beyond the RPA. It is well known that the ring ap- 
proximation can be applied satisfactorily for very small 
r s (r s < 1). In Fig. 3, however, we note that even for 
large r s , the variational rings give successful results. 

Here, we examined the efficiency of the VPT by cal- 
culating the correlation energy of the Coulomb electron 
gas. Further corrections beyond variational rings and the 
extension to finite temperatures are also easily accessible 
within the present VPT and we expect good applications 
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IV. SUMMARY 



Variational ring 
Conventional ring 
Monte Carlo 



After briefly reviewing the VPT formulated recently 
for the many-particle systems, we have applied this 
method to the Coulomb electron gas. The ground state 
correlation energy of the Coulomb electron gas is calcu- 
lated with the ring approximation of the VPT. Improve- 
ment is achieved by the expansion with variational prop- 
agators. Especially, we note that even for large r s , the 
approximation produces quite successful results, as is ev- 
ident from the good agreement with that of the Green's 
function Monte Carlo calculation. 



FIG. 3: The ground state correlation energy per particle cal- 
culated with ring diagrams; the solid and open circles are 
results of the variational and conventional ring approxima- 
tions, respectively. The triangles represent the result of the 
Green's function Monte Carlo calculation. 



to any other systems. 
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